###############################################################################################
# Effect of CRSV on Interethnic Relations
# 
# Figure G.1
# 
###############################################################################################

par(mfrow=c(1,3), mar=c(3, 5, 1, 0), cex.axis=.8, cex.lab=.8, cex.main=.8)


plot(c((inter.congo.ind), 
       (inter.congo.dir)), 2:1, pch=c(19, 21), cex=1.2, ylim=c(0.5, 2.5), xlim=c(-.6, .6), main="DRC", axes=F, ylab="", xlab="Difference in SD")

#grid(lty=1, lwd=.5, col="lightgrey")
abline(v=0, lty=2, col="darkgrey")
segments(c(inter.congo.ind.lo, 
           inter.congo.dir.lo), 2:1,
         c(inter.congo.ind.up, 
           inter.congo.dir.up), 2:1)
axis(2, 2:1, labels=c("List \n  Experiment", "Direct Item"), col="white")
axis(1, col="white")
points(inter.congo.dir, 1, pch=19, cex=.8, col="white")

text(c(mean(inter.congo.ind), 
       mean(inter.congo.dir)), 2:1, round(c(mean(inter.congo.ind), 
                                            mean(inter.congo.dir)),2), pos=3, cex=.8)
box(col="darkgrey")
par(mar=c(3, 2.5, 1, 2.5))

plot(c((inter.liberia.ind), 
       (inter.liberia.dir)), 2:1, pch=c(19, 21), cex=1.2, ylim=c(0.5, 2.5), xlim=c(-.6, .6), main="Liberia", axes=F, ylab="", xlab="Difference in Probability")


#grid(lty=1, lwd=.5, col="lightgrey")
abline(v=0, lty=2, col="darkgrey")
segments(c(inter.liberia.ind.lo, 
           inter.liberia.dir.lo), 2:1,
         c(inter.liberia.ind.up, 
           inter.liberia.dir.up), 2:1)
#axis(2, 3:1, labels=c("List \n  Experiment", "List Experiment \n (Misreporters)", "Direct Item"), col="white")
axis(1, col="white")
points(inter.liberia.dir, 1, pch=19, cex=.8, col="white")

text(c(mean(inter.liberia.ind), 
       mean(inter.liberia.dir)), 2:1, round(c(mean(inter.liberia.ind),
                                              mean(inter.liberia.dir)),2), pos=3, cex=.8)
box(col="darkgrey")
par(mar=c(3, 0, 1, 5))

plot(c(mean(inter.sri.ind), 
       mean(inter.sri.dir)), 2:1, pch=c(19, 21), cex=1.2, ylim=c(0.5, 2.5), xlim=c(-.6, .6), main="Sri Lanka", axes=F, ylab="", xlab="Difference in SD")

#grid(lty=1, lwd=.5, col="lightgrey")
abline(v=0, lty=2, col="darkgrey")
segments(c(inter.sri.ind.lo, 
           inter.sri.dir.lo), 2:1,
         c(inter.sri.ind.up, 
           inter.sri.dir.up), 2:1)
#axis(2, 3:1, labels=c("List Experiment", "List Experiment \n (Misreporters)", "Direct Item"), col="white")
axis(1, col="white")
points(inter.sri.dir, 1, pch=19, cex=.8, col="white")

text(c(mean(inter.sri.ind), 
       mean(inter.sri.dir)), 2:1, round(c(mean(inter.sri.ind), 
                                          mean(inter.sri.dir)),2), pos=3, cex=.8)

box(col="darkgrey")
################################################################################################
